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I. Introduction. 


This proposal requested funding to measure the durations of gamma-ray bursts (GRB) 
in the 4B catalog as well as to study the structure of GRB time profiles returned by the 
Burst And Transient Source Experiment (BATSE) on board the Compton Gamma- Ray 
Observatory. The duration (T90) was to be measured using the same techniques and 
algorithms developed by the principal investigator for the 3B data. The profile structure 
studies fall into the two categories of variability and fractal analyses. 


IIA. Duration Measurements 

As part of this contract, we agreed to measure the burst durations for the 4B catalog. The 
measurement of gamma-ray burst durations is significant for several reasons but especially 
so since duration is the only GRB characteristic which clearly shows that these objects fall 
into 2 distinct groups. Burst hardness also could be argued to be a parameter of equal 
importance, but a histogram of hardness alone does not show any grouping or clustering. 
However, we recently showed (Dezalay, Lestrade, et al, 1996: attached herein) that the 
hardness-duration diagram contains more information than previously believed. It was 
common knowledge that shorter bursts were harder, on average, than longer bursts. This 
anti-correlation has been well-documented. We found, however, that there is a secondary 
relationship hidden in the longer bursts. That is, we present clear evidence for a positive 
correlation in bursts whose durations (T90) are longer than 2 secs. This new finding has 
been confirmed by the work of Horack et al (1996). These new results make an accurate 
and dependable determination of T90 all the more important. 



IIB. Results of Duration Measurements 

There were a total of 429 bursts used in this study. These included trigger numbers between 
3177 and 5483 (inclusive). As was true in previous burst catalogs, data gaps (caused by 
the failure of the onboard tape recorders as well as other normal interruption of burst data 
accumulation) restricted the availability of some of the BATSE data types. Table I lists all 
429 bursts with the data types that were used to determine their durations. Appropriate 
comments will be made in the comments file, where the gaps restricted the T90 measure. 
The data type codes are l=discla, p=preburst, s=discsc, t=tte, and b=disclb. Disclb was 
used only to fill in gaps when discla data were unavailable. This was necessary for 65 
GRB’s. 

Two sample pages of output from our duration measuring program are included in this 
report as Figures 1-2. 

Of the 429. approximately 65 were short enough to require TTE data for accurate T90 
determination but their TTE datasets were missing packets. 

Figure 3 presents the duration histogram for these 429 GRB’s. As expected, the T90 
bimodality is still present. 

III. Fractal Algorithms 

Our current work in this interesting new area of mathematical physics concerns the ap- 
plication of several proven algorithms to GRB time profiles. In this effort I have the help 
of Van Yuan, a PhD graduate student at MSU (See attached resume). As part of his 
dissertation. Yan is looking for evidence of scale invariance embedded in the burst time 
profiles. The problem with studies of these kind is that noise and short data sets (both 
characteristics of GRB time profiles) obfuscate the results. However, that doesn t mean 
that we should give up. There are many different fractal algorithms - with different areas 
of applicability. We are searching for the correct algorithm to unlock the secrets of GRB. 
This search is especially important since standard methods of analysis have failed at almost 
every turn when applied to these cosmic events. 
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Table I: Data Types used for each Burst Profile. 

Trig# DLA Preb DSC TTE DLB Trig# DLA Preb DSC TTE DLB 


3177 

1 

P 

s 

- 

- 

3764 

1 

P 

s 

- 

- 

3178 

1 

P 

s 

- 

- 

3765 

1 

P 

s 

- 

- 

3193 

1 

P 

s 

- 

- 

3766 

1 

P 

s 

- 

- 

3200 

- 

P 

- 

- 

b 

3767 

- 

P 

s 

- 

b 

3212 

1 

P 

s 

- 

- 

3768 

1 

P 

s 

- 

- 

3215 

1 

P 

s 

- 

- 

3770 

1 

P 

s 

t 

- 

3217 

1 

P 

s 

- 

- 

3771 

1 

P 

s 

- 

- 

3218 

1 

P 

s 

t 

- 

3772 

1 

P 

s 

- 

- 

3220 

1 

P 

s 

- 

- 

3773 

1 

P 

s 

- 

- 

3227 

1 

P 

s 

- 

- 

3774 

1 

P 

s 

t 

- 

3229 

1 

P 

s 

- 

- 

3775 

- 

P 

- 

- 

b 

3237 

1 

P 

s 

- 

- 

3776 

1 

P 

s 

- 

- 

3238 

1 

P 

s 

- 

- 

3779 

1 

P 

s 

- 

- 

3240 

- 

P 

s 

t 

b 

3781 

1 

P 

s 

- 

- 

3241 

1 

P 

s 

- 

- 

3782 

1 

P 

s 

t 

- 

3242 

1 

P 

s 

- 

- 

3788 

1 

P 

s 

- 

- 

3243 

- 

P 

s 

t 

b 

3789 

1 

P 

s 

- 

- 

3245 

1 

P 

s 

- 

• 

3790 

- 

P 

s 

t 

b 

3246 

1 

P 

s 

- 

- 

3791 

1 

P 

s 

t 

- 

3247 

1 

P 

s 

- 

- 

3792 

1 

P 

s 

- 

- 

3248 

1 

P 

s 

- 

- 

3796 

- 

P 

- 

- 

b 

3249 

- 

P 

s 


b 

3797 

- 

P 

- 

t 

b 

3250 

1 

P 

s 

- 

- 

3799 

1 

P 

s 

t 

- 

3251 

1 

P 

s 

- 

- 

3800 

1 

P 

s 

- 

- 

3253 

1 

P 

s 

- 

- 

3801 

1 

P 

s 

- 

- 

3255 

1 

P 

s 

- 

- 

3803 

1 

P 

s 

t 

- 

3256 

1 

P 

s 

- 

- 

3804 

- 

P 

s 

t 

b 

3257 

1 

P 

s 

- 

- 

3805 

1 

P 

s 

- 

- 

3259 

1 

P 

s 

- 

- 

3806 

1 

- 

- 

- 

- 

3266 

1 

P 

s 

t 

- 

3807 

1 

P 

s 

- 

- 

3267 

1 

P 

s 

- 

- 

3810 

1 

P 

s 

t 

- 

3269 

1 

P 

s 

- 

- 

3811 

1 

P 

s 

- 

- 

3273 

1 

P 

s 

- 

- 

3812 

1 

P 

s 

- 

- 

3276 

1 

P 

s 

- 

- 

3814 

1 

P 

s 

- 

- 

3278 

1 

P 

s 

- 

- 

3815 

1 

P 

s 

- 

- 



Table I (cont.): Data Types used for each Burst Profile. 


Trig# 

DLA 

Preb 

DSC 

TTE 

DLB 

Trig# 

DLA 

Preb 

DSC 

TTE 

DLB 

3279 

1 

P 

s 

- 

- 

3819 

1 

P 

s 

- 

- 

3280 

1 

P 

s 

- 

- 

3840 

1 

P 

s 

t 

- 

3282 

1 

P 

s 

t 

- 

3843 

1 

P 

s 

- 

- 

3283 

1 

P 

s 

- 

- 

3846 

- 

P 

s 

t 

b 

3284 

1 

P 

s 

- 

- 

3848 

- 

P 

s 

t 

b 

3286 

1 

P 

s 

- 

- 

3853 

1 

P 

s 

- 

- 

3287 

1 

P 

s 

- 

- 

3857 

- 

P 

- 

- 

b 

3289 

- 

P 

s 

- 

b 

3860 

1 

P 

s 

- 

- 

3290 

1 

P 

s 

- 

- 

3864 

i 

P 

s 

- 

- 

3291 

- 

P 

- 

t 

b 

3866 

1 

P 

s 

- 

- 

3292 

1 

P 

s 

- 

- 

3867 

1 

P 

s 

t 

- 

3293 

1 

P 

s 

- 

- 

3868 

1 

P 

s 

- 

- 

3294 

1 

P 

s 

t 

- 

3869 

1 

P 

s 

- 

- 

3295 

- 

P 

s 

t 

b 

3870 

1 

P 

s 

- 

- 

3296 

- 

P 

s 

t 

b 

3871 

1 

P 

s 

- 

- 

3297 

1 

P 

s 

t 

- 

3875 

1 

P 

s 

- 

- 

3298 

- 

P 

- 

t 

b 

3879 

1 

P 

s 

- 

- 

3299 

- 

P 

s 

t 

b 

3885 

- 

P 

s 

- 

b 

3301 

1 

P 

s 

- 

- 

3886 

1 

P 

s 

- 

- 

3303 

1 

P 

s 

- 

- 

3887 

1 

P 

s 

- 

- 

3305 

- 

P 

s 

t 

b 

3888 

1 

P 

s 

t 

- 

3306 

1 

P 

s 

- 

- 

3889 

1 

P 

s 

t 

- 

3307 

1 

P 

s 

- 

- 

3890 

1 

P 

s 

- 

- 

3308 

1 

P 

s 

- 

- 

3891 

1 

P 

s 

- 

- 

3311 

- 

P 

- 

t 

b 

3892 

1 

P 

s 

- 

- 

3319 

1 

P 

s 

- 

- 

3893 

1 

P 

s 

- 

- 

3320 

1 

P 

s 

- 

- 

3894 

1 

P 

s 

t 

- 

3321 

1 

P 

s 

- 

- 

3895 

1 

P 

s 

- 

- 

3322 

1 

P 

s 

- 

- 

3897 

- 

P 

s 

- 

b 

3323 

1 

P 

s 

t 

- 

3899 

1 

P 

s 

- 

- 

3324 

1 

P 

s 

- 

- 

3900 

1 

P 

s 

- 

- 

3325 

- 

P 

s 

- 

b 

3901 

1 

P 

s 

- 

- 

3328 

- 

P 

s 

- 

b 

3902 

1 

P 

s 

- 

- 

3330 

1 

P 

s 

- 

- 

3903 

1 

P 

s 

- 

- 

3335 

1 

P 

s 

t 

- 

3904 

1 

P 

s 

t 

- 



Table I (cont.): Data Types used for each Burst Profile. 


Trig# DLA Preb DSC TTE DLB 

3336 1 p s 

3337 - p s - b 

3338 1 p s t 

3339 1 p s 

3340 1 p s t - 

3342 1 p s t 

3345 1 p s 

3347 1 p s - - 

3349 1 p s 

3350 1 p s - - 

3351 1 p s 

3352 1 p s 

3356 1 p s - - 

3357 1 p s 

3358 1 p s 

3359 1 p s t 

3360 p s - b 

3361 p s - b 

3364 1 p s 

3366 - p s - b 

3369 1 p s - - 

3370 1 p s 

3374 1 p s t 

3378 1 p s t 

3379 1 p s t 

3384 1 p s t 

3385 P s - b 

3403 1 p s 

3405 1 p s - - 

3406 1 p s 

3407 1 p s 

3408 1 p s 

3410 1 p s t 

3412 1 p s t 

3415 1 p s 


Trig# DLA Preb DSC TTE DLB 

3905 1 p s 

3906 1 p s - - 

3908 1 p s - - 

3909 1 p s 

3910 1 p s t 

3911 1 p s 

3912 1 p s 

3913 1 p s - - 

3914 1 p s 

3915 1 p s t - 

3916 1 p s - - 

3917 1 p s 

3918 1 p s 

3919 1 p s t 

3921 1 p s t - 

3924 1 p s 

3926 1 p s - 

3927 - p - - b 

3929 1 p s 

3930 1 p s 

3931 1 p s - 

3932 p - - b 

3935 1 p s 

3936 1 p s t 

3937 1 p s 

3938 1 p s 

3939 1 P s t - 

3940 1 p s - - 

3941 1 p s 

3954 1 p s - - 

4039 1 p s 

4048 1 p s 

4095 1 p s - 

4146 1 p s 

4157 1 p s - - 



Table I (cont.): Data Types used for each Burst Profile. 


Trig# 

DLA 

Preb 

DSC 

TTE 

DLB 

Trig# 

DLA 

Preb 

DSC 

TTE 

DLB 

3416 

1 

P 

s 

- 

- 

4216 

1 

P 

s 

- 

- 

3431 

1 

P 

s 

t 

- 

4251 

1 

P 

s 

- 

- 

3436 

1 

P 

s 

- 

- 

4256 

1 

P 

s 

- 

- 

3437 

1 

P 

s 

t 

- 

4312 

1 

P 

s 

- 

- 

3439 

1 

P 

s 

- 

- 

4327 

1 

P 

s 

- 

- 

3440 

- 

P 

s 

- 

b 

4350 

1 

P 

s 

- 

- 

3441 

1 

P 

s 

t 

- 

4368 

1 

P 

s 

- 

- 

3442 

- 

P 

- 

- 

b 

4388 

1 

P 

s 

- 

- 

3443 

- 

P 

- 

- 

b 

4462 

- 

P 

s 

- 

b 

3443 

1 

P 

s 

- 

- 

4469 

1 

- 

- 

- 

* 

3449 

- 

P 

s 

t 

b 

4556 

1 

P 

s 

- 

- 

3450 

- 

P 

- 

t 

b 

4569 

1 

P 

s 

- 

- 

3453 

1 

P 

s 

- 

- 

4636 

- 

P 

s 

- 

b 

3458 

1 

P 

s 

- 

- 

4649 

- 

P 

- 

- 

b 

3464 

1 

P 

s 

- 

- 

4653 

1 

P 

s 

- 

- 

3465 

1 

P 

s 

- 

- 

4660 

1 

P 

s 

t 

- 

3466 

1 

P 

s 

- 

- 

4701 

1 

P 

s 

- 

- 

3467 

1 

P 

s 

- 

- 

4710 

1 

P 

s 

- 

- 

3471 

1 

P 

s 

- 

- 

4744 

1 

P 

s 

- 

- 

3472 

1 

P 

s 

- 

- 

4745 

1 

P 

s 

- 

- 

3473 

1 

- 

s 

- 

- 

4757 

- 

P 

- 

- 

b 

3476 

1 

P 

s 

- 

- 

4761 

- 

P 

- 

t 

b 

3477 

1 

P 

s 

t 

- 

4776 

1 

P 

s 

- 

- 

34S0 

1 

P 

s 

t 

- 

4807 

1 

P 

s 

- 

- 

3481 

1 

P 

s 

- 

- 

4814 

1 

P 

s 

- 

- 

3485 

1 

P 

s 

- 

- 

4871 

1 

P 

s 

- 

- 

3486 

1 

P 

5 

- 

- 

4S9S 

- 

P 

- 

- 

b 

3487 

1 

P 

s 

t 

- 

4939 

1 

P 

s 

- 

- 

34SS 

1 

P 

s 

- 

- 

4955 

1 

P 

s 

t 

- 

3489 

1 

P 

s 

- 

- 

4959 

1 

P 

s 

- 

- 

3491 

1 

P 

s 

- 

- 

5079 

1 

P 

s 

t 

- 

3492 

1 

P 

s 

- 

- 

5080 

1 

P 

s 

- 

- 

3493 

1 

P 

s 

- 

- 

5123 

1 

P 

s 

t 

- 

3494 

1 

P 

s 

- 

- 

5206 

1 

P 

s 

t 

- 

3499 

- 

P 

s 

- 

b 

5212 

1 

P 

s 

t 

- 

3502 

1 

P 

s 

t 

- 

5255 

1 

P 

s 

- 

- 



Table I (cont.): Data Types used for each Burst Profile. 


Trig# 

DLA 

Preb 

DSC 

TTE 

DLB 

Trig# 

DLA 

Preb 

DSC 

TTE 

DLB 

3503 

1 

P 

s 

- 

- 

5277 

1 

P 

s 

t 

- 

3505 

1 

P 

s 

- 

- 

5299 

1 

P 

s 

- 

- 

3509 

1 

P 

s 

- 

- 

5304 

1 

P 

s 

- 

- 

3510 

1 

P 

s 

t 

- 

5305 

1 

P 

s 

- 

- 

3511 

1 

P 

s 

- 

- 

5316 

1 

P 

s 

t 

- 

3512 

1 

P 

s 

- 

- 

5337 

1 

- 

- 

- 

- 

3514 

1 

P 

s 

- 

- 

5339 

1 

P 

s 

- 

- 

3515 

1 

P 

s 

- 

- 

5377 

- 

P 

- 

- 

b 

3516 

1 

P 

s 

- 

- 

5379 

1 

P 

s 

- 

- 

3523 

1 

P 

s 

- 

- 

5387 

1 

P 

s 

- 

- 

3527 

1 

P 

s 

- 

- 

5389 

1 

P 

s 

- 

- 

3528 

1 

P 

s 

- 

- 

5392 

- 

P 

* 

- 

b 

3529 

- 

P 

s 

- 

b 

5407 

1 

P 

s 

- 

- 

3530 

1 

P 

s 

- 

- 

5409 

1 

P 

s 

- 

- 

3537 

- 

P 

s 

- 

b 

5410 

1 

P 

s 

- 

- 

3545 

1 

P 

s 

t 

- 

5411 

1 

P 

s 

- 

- 

3552 

1 

P 

s 

- 

- 

5412 

1 

P 

s 

- 

- 

3567 

1 

P 

s 

- 

- 

5413 

1 

P 

s 

- 

- 

3569 

1 

P 

s 

- 

- 

5415 

1 

P 

s 

- 

- 

3571 

1 

P 

s 

- 

- 

5416 

1 

P 

s 

- 

- 

3580 

1 

P 

s 

- 

- 

5417 

1 

P 

s 

- 

- 

3585 

1 

P 

- 

t 

- 

5418 

- 

P 

s 

- 

b 

3588 

1 

P 

s 

- 

- 

5419 

1 

P 

s 

- 

- 

3590 

- 

P 

- 

- 

b 

5420 

1 

P 

s 

- 

- 

3593 

1 

P 

s 

- 

- 

5421 

1 

P 

s 

- 

- 

3594 

1 

- 

s 

- 

- 

5423 

1 

P 

s 

- 

- 

3595 

- 

P 

- 

- 

b 

5425 

1 

P 

s 

- 

- 

3598 

1 

P 

s 

- 

- 

5427 

1 

P 

s 

- 

- 

3606 

1 

P 

s 

t 

- 

5428 

1 

P 

s 

- 

- 

3608 

1 

P 

s 

- 

- 

5429 

1 

P 

s 

- 

- 

3611 

i 

i 

P 

s 

- 

- 

5432 

- 

P 

- 

- 

b 

3618 

1 

P 

s 

- 

- 

5433 

1 

P 

s 

- 

- 

3634 

1 

P 

s 

- 

- 

5434 

1 

P 

s 

- 

- 

3637 

1 

P 

s 

- 

- 

5436 

1 

P 

s 

- 

- 

3639 

1 

P 

s 

- 

- 

5439 

1 

P 

s 

- 

- 

3640 

1 

P 

s 

t 

- 

5443 

- 

P 

s 

- 

b 




Table I (cont.): 

Data Types used for each Burst Profile. 


Trig# 

DLA 

Preb 

DSC 

TTE 

DLB 

Trig# DLA Preb 

DSC 

TTE 

DLB 

3642 

1 

P 

s 

- 

- 

5444 1 

P 

s 

- 

- 

3643 

1 

P 

s 

t 

- 

5446 1 

P 

s 

- 

- 

3644 

1 

P 

s 

t 

- 

5447 1 

P 

s 

- 

- 

3647 

1 

P 

s 

- 

- 

5448 1 

P 

s 

t 

- 

3648 

1 

P 

s 

- 

- 

5450 1 

P 

s 

- 

- 

3649 

1 

P 

s 

- 

- 

5451 1 

P 

s 

- 

- 

3651 

- 

P 

- 

- 

b 

5452 1 

P 

s 

- 

- 

3652 

1 

P 

s 

- 

- 

5453 1 

P 

s 

t 

- 

3654 

1 

P 

s 

- 

- 

5454 1 

P 

s 

- 

- 

3655 

1 

P 

s 

- 

- 

5456 1 

P 

s 

- 

- 

3657 

- 

P 

- 

- 

b 

5457 1 

P 

s 

- 

- 

3658 

1 

P 

s 

- 

- 

5458 1 

P 

s 

t 

- 

3662 

1 

P 

s 

- 

- 

5459 1 

P 

s 

t 

- 

3663 

1 

P 

s 

- 

- 

5461 1 

P 

s 

- 

- 

3664 

1 

P 

s 

- 

- 

5462 1 

P 

s 

- 

- 

3665 

1 

P 

s 

t 

- 

5463 1 

P 

s 

- 

- 

3668 

1 

P 

s 

t 

- 

5464 1 

P 

s 

- 

* 

3671 

1 

P 

s 

- 

- 

5465 1 

P 

s 

t 
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ABSTRACT 

The hardness-duration diagram for gamma-ray bursts still contains undiscovered important information about 
the intrinsic properties of these phenomena. We analyze these diagrams for the PHEBUS and BATSE 
experiments First, we show that the BATSE diagram is very similar to that for PHEBUS when we restrict the 
BATSE data set to events observable by PHEBUS. In this case, both diagrams present a high degree of clustering 
into two subclasses. This shows that the brightness of the events is a more important factor in determining the 
aspect of this diagram than the hardness ratio energy ranges. Second, for the subclass of long bright bursts, both 
experiments show evidence for a positive correlation between hardness and duration. This is a significant new 
result, as it represents an intrinsic property of long events. The commonly held perception of an anticorrelation 
(between hardness and duration) for all bursts must be replaced by a more complicated model that treats the two 
subclasses separately. No statistically significant correlation is found within the short-event population. The 
existence of such an intrinsic positive correlation for the long bursts would complicate the search for cosmologica 
effects in the gamma-ray burst data. 

Subject heading: gamma rays: bursts 


1. INTRODUCTION 

The most significant advances in our understanding of the 
nature of cosmic gamma-ray bursts (GRBs) have been made 
through statistical studies. Statistical studies of the position 
and of the intensity of bursts have restricted models of the 
spatial distribution of sources to either an extended galactic 
halo or a cosmological distribution. They are also the principal 
tool used to search for subclasses in the GRB population. To 
date, morphological studies of burst time profiles have not 
succeeded in exhibiting any classification schemes (Fishman & 
Meegan 1995; Lestrade 1994). The only strong evidence that 
subclasses of events exist comes from the durations of GRBs. 
The first piece of evidence is the bimodality of the duration 
distribution reported by several authors (Klebesadel 1992; 
Kouveliotou et ai. 1993). The second is that the spectrum of 
short events is, on average, harder than that for longer events 
(Dezalay et al. 1992, 1995, 1996; Kouveliotou et al. 1993). At 
present no significant differences have been detected between 
the angular or intensity distributions of the two kinds of bursts. 
The only evidence we have that the population of GRBs may 
be divided into at least two subclasses lies in the hardness- 
duration diagram (HDD). This diagram could, therefore, be 
considered as a potential Hertzspmng-Russeii diagram for 
GRBs. In the following, we compare the HDD from BATSE 
and PHEBUS, which are the only experiments presenting a 
significant difference between the average hardness of short 
and long bursts. We perform this analysis using the Third 
BATSE Catalog (Meegan et al. 1996) and the PHEBUS data 
set. We also study the hardness and duration distributions of 
short-anri Jong evens se pai aic iy in order to seared ior spcc rik 
trends. 


2. THE HARDNESS-DURATION DIAGRAMS 

The qualitative comparison of the PHEBUS and BATSE 
HDDs gives the impression that the short-burst and long-burst 
populations are more separated in PHEBUS than they are in 
BATSE data. Although not statistically significant, this fact has 
to be investigated further, since one would have expected a 
larger separation in BATSE due to the stronger bimodalicy of 
the duration distribution. It is important to determine whether 
the different aspects of the two diagrams are due to the energy 
bands used to calculate the hardness, to the intensity of events, 
or to a smaller event sample in PHEBUS. 

We calculated the durations of PHEBUS (Barat et al, 1988) 
events using the algorithm by Kouveliotou et al. (1993). 
The major advantage of this integral criterion is that it gives 
measures that do not depend on the background level. Other 
measures were tried, including a differential one based on the 
intensity of the time profile above the background level 
(Dezalay et al. 1995). We find that the choice of the duration 
criterion has no effect on the shape of the the HDD or on the 
conclusions drawn from it. The bimodality is only slightly 
present in the distribution of PHEBUS 7*} durations. This fact 
shows that the separation observed for the two HDD popula- 
tions in PHEBUS is not the consequence of a strong bimodai- 
ity. The hardness of PHEBUS events is calculated using the 
mean hardness ratio (MHR). It is derived from the count 
spectrum integrated over the 3 cr duration defined in Dezalay 
et al. (1995). The energy ranges used are 120-320 keV and 
320-7000 keV. The distribution of the mean hardness ratios is 
unimodal. The HDD constructed using MHR and 7** is 
reported in Dezalav et al. (1996) and displayed in Figure la. 

To compare PHEBUS and BATSE, we use the data pub- 
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178 type l events are not displayed. 

lished in the Third BATSE Gamma-Ray Burst Catalog (3B 
Catalog) by Meegan et al. (1996). The hardness of BATSE 
events is measured using the ratio of fluences from different 
channels. Uncertainties associated with this quantity are larger 
than those of the hardness ratio (HR) of counts. To illustrate 
the difference between the two criteria, we compare the 
significance of the separation between the average hardness of 
short (T90 < 2 s) and long (7, 0 > 2 s) events of the IB 
Catalog. Using HR 32 (ratio of the total counts in the 100- 
300 keV energy range over those in the 50-100 keV range), 
Kouveliotou et al. (1993) reported that the difference in 
hardness is significant at a 7.3 cr confidence level. The same 
quantity calculated with the ratio of fluence 3 (100-300 keV) 
to fluence 2 (50-100 keV), ER, :> is only 2.2 cr. This result is 
explained mostly by the large uncertainties of the fluences of 
short bursts due to the deconvolution process on a small 
number of counts. The relative uncertainty of the average 
hardness of short events is 5% using HR 32 , while it is 25% 
when calculated with FR 32 . 

In order to put the two experiments on an equal footing of 
sensitivity, we have to select only the most intense BATSE 
events. To this end, we used the peak flux on the 256 ms 
timescale (Piss) in the 50-300 keV energy range as our 
intensity criterion for BATSE. This timescaie corresponds to 
the shortest trigger time interval in PHEBUS. Selecting the 
most intense bursts for BATSE according to Fist will affect 
events shorter than j s the same way in BATSE as in 
PHEBUS. The 3B Catalog and the PHEBUS data set have 51 
events m common. Among this sample, only 39 have an 
available measure of three quantities: T w , fluence, and peak 


flux. The weakest of these 39 events has a peak flux equal to 
2.8 photons cm" 1 s' 1 . In this paper we refer to the sample of 
BATSE bright bursts with a peak flux above this value as type 
I. It contains 178 events. Similarly, we construct two other 
BATSE samples according to their intensity: type II comprises 
446 events with 0.6 < Pus < 2.8 photons cm' 2 s' 1 , and type 
HI has 178 events with Pus <0.6 photons cm' 2 s" 1 . 

As presented in Figure 1 b, the HUD of BATSE type I 
bursts is qualitatively very similar to the PHEBUS HDD (Fig. 
la). To be able to compare the diagrams by eye more easily, 
we have smoothed the imaginary line linkin g all points (sorted 
according to their duration) using a weighted Loess filtering 
method. The trend of this smoothed line is similar for the 
BATSE type I and PHEBUS bursts and consists qualitatively 
of two connected power laws, the first with a negative power- 
law index up to T*, ~ 5 s, followed by the second with a 
positive slope. Although no quantitative information can be 
extracted from this smoothed line, it shows that the diagrams 
look similar even when different criteria are used to measure 
the hardness and different energy ranges to measure the 
duration. The positive trend observed within the long bursts 
suggests that a separate analysis for the two populations might 
prove fruitful. 

3. CORRELATION BETWEEN HARDNESS AND DURATION 

In this section we restrict our analysis to samples of long 
events (7« 0 > 2 s) detected by BATSE and PHEBUS. We 
calculate the significance of the positive trend observed in the 
HDD using two different statistical tests. We investigate the 
evolution of this correlation between spectral hardness and 
duration as a function of the intensity of bursts in the sample 
in BATSE. We also briefly discuss systematic effects that could 
affect these observations. 

3.1. Positive Correlation Hardness Ratio and Too in 
Long and Bright Bursts 

Long and bright bursts (7, 0 >2 s) from the BATSE type I 
and PHEBUS samples are displayed in the HDD in Figure 2. 
The significance of the trend observed in Figure 2 is estimated 
by the Spearman and the Kolmogorov-Smimov (K-S) tests 
(Press et al. 1992). The Spearman rank-order test applied to 
PHEBUS long bursts (see Table 1) yields a correlation 
coefficient (Spearman rho) equal to 0.243. The probability 
(•Pspcr) that this correlation is due to a random fluctuation is 
5.1 X 10 ~ 3 and corresponds to a 2.8 cr confidence level. The 
same test applied to the BATSE type I events yields a similar 
result with a correlation coefficient of 0.216, or 2.6 cr 
t P . = 9.7 x 10" 3 ). Since the two results are independent, the 
probability of finding the same correlation in both PHEBUS 
and BATSE is equal to 5 x 10 5 . 

Dividing the long bursts into two subsamples according to 
the median duration allows us to apply the K-S test. The 
hypothesis tested is that the hardness distributions of the two 
subsamples come from the same parent population. The K-S 
test confirms the result obtained with the Spearman test, i.e., 
that in bright and long bursts the spectral hardness is positively 
correlated with the duration. Results of both tests applied to 
the long bursts in each of the three BATSE samples and 
PHEBUS are summarized in Table 1. 

If there is a correlation between 7*, and brightness (Norris 
et al. 1996) associated with the well-known correlation be- 
tween hardness and brightness (Atteia et al. 1994, Nemiroff et 
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T 90 (s) 

Fig. 2. — HDD for long (r w >2 s) PHEBUS bursts ( filled circles) and 
BATSE type I bursts (open circles). The two lines represent the least-squares 
fits for BATSE (solid line ) and PHEBUS ( doited line) events in the HDD, 


Peak flux (ph.cm^.s' 1 ) 

Fig. 3. — Probability P Spe „ of obtaining, by chance, a positive correlation 
between FR^ and Too for a sample of lone BATSE bursts brighter than the 
threshold value of the peak flux (abscissa) on the 256 ms ( solid line) and 
1024 ms ( dotted line) timescales. 


al. 1994), it could result in a positive correlation between 
duration and hardness. However, a recent study bv Koshut et 
a I (1956) of systematic effects on 7^ shows that there is no 
systematic dependence of 7 90 on the burst signai-to-noise ratio 
for a sample of simulated and real time profiles. Moreover, if 
there is a correlation between 7^ and intensity, it should also 
be observed in weaker burst samples when the brightness 
range is the same. This is not the case (see § 3.2). Therefore, 
since the correlation is still observable and even more signif- 
icant using the brightest BATSE events, we think this result 
cannot be explained by selection or systematic effects in the 
data. We believe that we have identified for the first time an 
intrinsic property within the long-burst population. We are not 
aware of any theoretical model predicting a positive correla- 
tion between hardness and duration. On the other end, some 
models, e.g., Meszaros & Rees (1994), predict an anticorreia- 
tion between hardness and duration. However, there is no 

TABLE 1 

Statistics of the Four Samples 


Statistic PHEBUS Type l Type II Type III 

Number short/long 1 * 43/131 34/144 139/307 35/143 

Separation in hardness b 6.0 4.7 5.5 13 

Spearman rmr 0.243 0.216 -0.109 0.063 

P Speaf 0.0051 0.0097 0.056 0.46 

/> w 0.046 0.0060 0.088 0.30 

With equal brightness range in BATSE 

Number short/long* 9/44 57/179 142/371 

Spearman rho c *** 0355 0.020 —0.027 

p Spat 9.1 x 10" 5 0.79 0.60 


1 Using the Too 3 * 2 s boundary. 

^iumfieance m me sep arat ion o e t wegm fee avera ge hanmeSROtfiTTOTI <<2 s) 

and long (>2 s) events. 

c Spearman rank-order correlation coefficient. 


contradiction since they do not treat the two subclasses 
separately. 

3.2. Evolution of the Correlation With Intensity 

We do not contend that the whole ensemble of bursts 
displays a positive correlation between burst hardness and 
duration. On the contrary, the correlation we find exists for 
only the nearest and brightest long GRBs. No significant 
correlation has been found in the weaker type II and III events 
(Table 1), or in the short-burst population. The three samples 
of bursts in BATSE were created to compare PHEBUS and 
BATSE. However, the range in the intensities of type I events 
is larger than that in type II and III samples. To check whether 
different ranges could affect the correlation, we performed the 
same analysis with three new samples having comparable 
ranges in brightness. Results summarized in Table 1 confirm 
that the correlation is still not observed in the two weaker 
event samples. 

To take this analysis one step further, we have calculated the 
probability of finding a correlation in the long BATSE bursts 
by chance using different peak flux criteria. Figure 3 shows that 
as the peak flux threshold (on the 1024 and 256 ms timescales) 
is increased (meaning that we are selecting brighter and 
brighter bursts), the correlation between hardness and dura- 
tion improves. The lowest probability observed is 10 ~ 5 . It is 
calculated using the ~50 brightest long BATSE events. When 
the value of ~7 photons cm * s" 1 is reached for peak flux, only 
46 long events remain. Past this value, the increase in proba- 
bility is likely due to poor statistics. We have repeated this 
calculation for the peak flux integrated on 64 ms and found the 
same behavior. Figure 3 generalizes the results shown in Table 
1, that the correlation increases with the average brightness of 
the long bursts. This evolution could be due to cosmological 
effects becoming more important for weaker events. At this 
time, we cannot exclude the possibility that the weakening of 
the correlation between spectral hardness and duration ob- 
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served for, long and bright bursts is caused by systematic 
effects. 

4. CONCLUSIONS 

We have shown that the HDD diagrams from PHEBUS and 
from a subsampie of BATSE bright events are similar. This 
result tends to show that the difference in the energy ranges 
used to measure the hardness between the two experiments 
has very little influence on the aspect of the diagram. We have 
also found a positive correlation between the duration and the 
hardness in the samples of long and bright bursts detected by 
PHEBUS and BATSE. This correlation is not detectable in 
samples of long and weak events or in the short-event popu- 
lation. The different trends exhibited by the two subclasses in 
the diagram suggest that they should be studied separately. 
This result has an important implication for the research of 
cosmological effects in the data. For an extragalactic distribu- 
tion of sources, we expect longer bursts to be weaker due to 


time dilation, and weaker events to be softer due to the 
redshift of the spectrum (Paczynski 1992). Therefore, since we 
find that within the long-burst population the longer events are 
harder on average, this result cannot be explained by cosmo- 
logical effects. We believe that the correlation is an intrinsic 
property of the long-burst population. The presence in the 
data of such a correlation does not contradict or support the 
cosmological hypothesis, since it is observable only for bright 
events. Nevertheless, this correlation between hardness and 
duration may complicate the search for cosmological effects. 

The positive correlation (hardness vs. duration) has been 
confirmed independently by Horack et al. (1996) using the 50 
brightest bursts in the BATSE 3B Catalog. 

We thank Tom Koshut (USRA) for fruitful discussions. We 
also acknowledge the help of James Berger (Purdue Univer- 
sity) for his clustering analysis. J. P. L. was partially supported 
by NASA grant NCC8-82. 
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